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Abstract: It is shown that the " chaoticity hypothesis" , analogous to Ruelle's principle for tur- 
bulence and recently introduced in statistical mechanics, implies the Onsager reciprocity and the 
fluctuation dissipation theorem in various reversible models for coexisting transport phenomena. 



§1: Introduction. 

In reference [GC2] we introduced, as a principle holding when a system motions have an 
empirically chaotic nature, that: 

Chaotic hypothesis: A many particle system in a stationary state can be regarded as a smooth 
dynamical system with a transitive axiom A global attractor for the purpose of computing macro- 
scopic properties. In the reversible case it can be regarded, for the same purposes, as a smooth 
transitive Anosov system. 

See [AA],[S],[Rl],[Bo] for the notion and properties of Anosov systems and [Bo,R2] for the more 
general systems with Axiom A attractors. I do not know examples of reversible systems with 
an Axiom A transitive global attractor which are not transitive Anosov systems: in this sense 
I regard the second part of the hypothesis as a likely consequence of the first (the viceversa is 
trivial). In this paper only the latter part of the hypothesis will be needed and used. 
The flow x —>■ Vtx solving in phase space the differential equations of motion generates an 
evolution t — ► F(Vtx) on the observables F(x). It can be naturally interpreted as a random 
process when the initial dat chosen randomly with a given distribution /io- 

The averages over the time variable t give the stationary state for the evolution Vt, provided 
they are uniquely defined, i.e. provided for almost all choices of x with distribution /io the 
averages exist and are ^-independent. In this case the stationary state is a stationary probability 
distribution [i: to stress that it is dependent on /io we ca U it the statistics /j, of [lo- 
in many applications it is more convenient to regard the evolution as a discrete trasformation 
defined on a restricted phase space C of observed events, also called timing events, (which could 
be, for instance, the occurrence of a microscopic binary "collision"). The time evolution, or 
the dynamics, is a map S of C into itself. The map S is derived from the flow Vt solving 
the differential equations of motion of the system: the timing events C have to be thought as 
a surface transversal to the flow and if t(x) is the time between the timing event x and the 
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successive one Sx it is V t ( x )X = Sx. Note that the points Vtx are not timing events (i.e. they 
are not in C) for the intermediate times < t < t(x). 

The notions of statistics /j, of po carry over unchanged to the discrete notions of phase space 
and of evolution and the above chaotic hypothesis is assumed in such a context. When referring 
to phase space, unless stated otherwise, we think of a phase space C consisting of timing events 
and of a map S on C defining the time evolution. The smoothness of the Anosov system is 
intended in all the coordinates and parameters on which the system equations depend. 

The chaotic hypothesis implies, as a mathematical consequence, that for most random distri- 
butions no for the choice of the initial data x the distribution /j, exists. 

However the choice of the initial data with distribution po proportional to the volume mea- 
sure on C plays a special role, because in the case of hamiltonian systems such distribution is 
generated naturally via the microcanonical ensemble. 

For instance one can read, translating symbols into the present notations: "the appropriate ob- 
jects of study of a statistical mechanics of time dependent phenomena are the random processes 
F(Vtx) with initial distribution of x, Pe(%) (the microcanonical distribution), for all energies of 
interest and for all gross variables F of interest" , (italics added), see the nice paper [GCDR]. 

I "strongly disagree" with such a philosophical position: but I will adopt it here because it 
is obvious that, whether one agrees or not, it is of fundamental interest to understand the 
statistics of initial data chosen at random with a distribution proportional to the phase space 
volume. Other random choices may have to wait until the latter is properly understood. 

In the Physics literature the existence of the distribution p is, in fact, assumed in general 
as stated by the following (extension) of the zeroth law, [UF], giving a global property of the 
motions generated by initial data chosen randomly with distribution po proportional to the 
volume measure on C: 

Extended zero-th law: A dynamical system (C,S) describing a many particle system (or a con- 
tinuum such as a fluid) generates motions that admit a statistics p in the sense that, given 
any (smooth) macroscopic observable F defined on the points x of the phase space C, the time 
average of F exists for all no -randomly -chosen initial data x and is given by: 

1 M_1 t 
J™, M £ F ^ x ) = / ^x')F{x') (1.1) 

k = Jc 

where p is a S-mvariant probability distribution on C. 

The chaotic hypothesis was proposed by Ruelle in the case of fluid turbulence, and it is extended 
to non equilibrium many particle systems in [GC1]. If one assumes it, then it follows that the 
zeroth law holds, [S,Bo,R]; however it is convenient to regard the two statements as distinct 
because the hypothesis we make is " only" that one can suppose that the system is Axiom A 
or Anosov for "practical purposes": this leaves the possibility that it is not strictly speaking 
such and some ("neglegible in the thermodynamic limit") corrections may be needed on the 
predictions obtained by using the hypothesis. 

From now on only reversible systems will be considered in this paper: they are dynamical 
systems such that there is an isometric map i of phase space such that i 2 = 1 and iS = S~ 1 i. 
In [GC2] the generality of the hypothesis is discussed and in [GC1], [GC2] we derived, as a 
rather general consequence, predictions testable by numerical experiment in systems with few 
degrees of freedom. The most relevant feature of the prediction, which is a large deviation 
theorem or fluctuation theorem for systems with reversible dynamics, is that it is parameter 
free. 

A drawback is that it is not testable directly in really large systems. 

The question of whether the "chaotic hypothesis" could be tested on real experiments, i.e. for 
really macroscopic systems, remained open. In this paper we show, through examples, that the 
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chaotic hypothesis implies quite generally, in systems with reversible dynamics, the Onsager 
reciprocity and the fluctuation-dissipation theorem: see [DGM] for a classical general discussion 
of the reciprocity relations, see [C] (Appendix A, p. 187-200) for a kinetic derivation, see also 
[E,ELS, GJL] for recent developments. A very nice introduction and an exposition of the basic 
ideas can be found in [D]. The ideas of the present paper can be applied also to systems relevant 
for the theory of developed turbulence (see [G4]). 
A puzzling aspect of the chaotic hypothesis is that it implies that the system has a positive 
gap separating from the Lyapunov exponents and one may have serious doubts on the validity 
of such a strong property: so that a discussion is in order. 

This is discussed in [GC2], §8, and [G4] by suggesting that there may well be many vanishing 
exponents (or exponents of order 0(N~ 1 )) if N is the particle number: such exponents should 
be "ignored" as they should correspond to macroscopic evolution laws which microscopically will 
be effectively described by local conservation laws [EMY,KV,LOY]. They have an approximate 
character, unless N —> oo, but one can think of imposing them as exact conservation rules by 
adding to the forces acting on the system other suitable "auxiliary" forces minimally required to 
achieve the purpose of turning the slow macroscopic observables responsible for the existence 
of the Lyapunov exponents into exact local conservation rules. For instance one can find 
the auxiliary forces by applying the Gauss' principle of least constraint, see appendix. The 
dynamical system obtained in this way should be one to which the chaotic hypothesis should 
apply. 

In §2, §3 we introduce some models which will be used to illustrate our general ideas; in §4 we 
discuss the relevant mathematical facts about reversible Anosov systems; in §5 we give a proof, 
whose full mathematical rigor still rests on a mathematical conjecture (§5) on Anosov systems 
(that I hope to address elsewhere), of the validity of the fluctuation dissipation relation and 
of Onsager reciprocity in the models introduced in §2, §3: but the generality of the argument, 
which seems largely model independent, will also emerge. Comments and a critical comparison 
with the classical derivations (that apply to our models as well) are presented in §6. 

A diffusion problem. 

We consider a mixture of two chemically inert gases whose N = N\ + N2 molecules (with 
equal masses m and respective numbers N\ = N2) are contained in a box B = [—^L, \L] 2 
with periodic boundary conditions and are subject to a respective force field E} = E 1 i and 
E_ 2 = E 2 i. The molecules interact via a pair interaction with a short range potential {e.g. a 
Lennard Jones type of potential). 

Furthermore the motion is subject to the constraint that the total energy is constant, via 
a constraint force law which is ideal in the sense that it verifies Gauss' principle of minimal 
constraint. This means that, if Ej is the field on the j-th particle (equal either to E 1 or to E 2 ) 
the equations of motion are: 

q . = — p ., p . = F ; + Ej i — ap . (2.1) 
—3 m—y —3 ~ 3 ~ 3 v ' 

with a = — ^ and we expect that the future time average of the total momentum P_ 

^3 —j 

and of the dissipation a will be some (P_)+ > and (a)+ > if (E 1 ^ 2 ) ^ 0. 
The average is expected to be attained exponentially fast while the "conjugate" variable C_, 
the center of mass position, will be expected to evolve with zero Lyapunov exponent (and to 
behave asymptotically as (C_)+ = (jVt+jv^m (H.)+t + cost). 

According to the analysis of [GC2],§8, and [G4], we expect that we can freely add to the 
equations of motion a minimal constraint force imposing the constraint P = P(E 1 ,E 2 )i if 
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P(E 1 ,E 2 ) is the average horizontal momentum, i.e. the (unknown) "equation of state", (and 
the implied C_ = cost). This means that considering a modified equation of motion: 

q. = —p., P ■ = F , + Eji - ap . - (3 (2.2) 

with /3 = ' El n\ Xn 3 2 Ni ^ ~ 01 n^n 2 ' s h°uld not lead to appreciably different qualitative behaviour 
if the initial data are consistent with the equation of state P_ = P(E 1 , E 2 ) i, C_ = 0. 
The phase space contraction per unit time is: 

1 



7 = (2(#! + N 2 )(l 



■)-l)a = (2(N 1 + N 2 )-3)a 



(2.3) 



#i +N 2 

Note that without imposing the extra constraint the phase spave contraction would have been 
(2(#i + N 2 ) — I) a: with a relative difference of 0(N~ 1 ), which should become neglegible in 
the thermodynamic limit L —> oo, N\L~ 2 = pi, N 2 L~ 2 = p 2 . 

We call Q the work per unit time performed by the forces ap^ so that the phase space con- 
traction rate in the configuration x is (to leading order in N) 2Na{x), i.e. E Niai ^l^ x ^ 2a2 ^ 

= jQjt if ks is Boltzmann's constant, ksT(x) is the kinetic energy per particle and ai(x), a 2 (x) 
are the horizontal coordinates of the centers of mass C_ 1 , C_ 2 of the two species, and ai(x), a 2 (x) 
the corresponding velocities. Thus the phase space contraction rate can be interpreted as 
times the entropy creation rate (at least if E 1 , E 2 are so small that the centers of mass horizontal 
velocities are small compared to the root mean square velocities). 

The above model is closely related to the color diffusion model considered in [BEC]. A very 
important feature of the model is its time reversibility: the map i(p, q) = (— p, q) has the 
property that it is an isometry of phase space such that iS = S~ 1 i and i 2 = 1. 

§5: A heat conduction-electrical conduction model. 

As a second model we consider a modification of model 4 of [CG2], inspired by [HHP], see 
also [PH] for a more general perspective on constrained systems. This is a model for a heat 
conducting and electrically conducting gas. In a box B = [— 2L — H ,2L + H] x [—^L, ^L] are 
enclosed N particles with mass m, interacting via a rather general pair potential, like a hard 
core potential with a tail or via a Lennard Jones potential, and they are subject to a constant 
force field (electric field) E i in the x direction. The boundary conditions are periodic in the 
horizontal direction and reflecting in the vertical direction. 

Adjacent to the box B there are two boxes 1Z+,1Z- containing N + = N_ particles interacting 
with each other via a hard core interaction, and with the particles in B via a pair interaction 
(e.g. with potential equal to the one between the particles in B), but are separated from the 
latter by a reflecting wall. 




The model name is motivated because we imagine other forces to act on the system: they 
are the minimal forces (in the sense of Gauss' minimal constraint principle, see appendix Al) 
necessary to enforce the following constraints: 
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1) the total kinetic energy in the "hot plate" 1Z+ and in the "cold plate" 1Z- are constrained 
to be N + k B T + and, respectively, N_k B T_ where T_ and T+ = T_ + 6T, 8T > are the 
temperatures of the plates. 

2) the total energy U in the box B is constrained to stay constant. 
The equations of motion are: 

P ■ 
—j 

-i ~ m (3.1) 

Pj = Fj + E x(qj)i ~ a+X+ig^Pj ~ a-X-ig^Pj ~ axig^Pj 

where x>X± are the characteristic functions of the regions B,1Z± and a + ,a_,a are multipliers 
defined so that for some T±,U: 

N± p 2 N p 2 

X>±(iP^ = N ± k B T ± , + V(q) = U (3.2) 

i=i i=i 

are exact constants of motion. 

We suppose for simplicity (see however comment (1) in §6 below) that the system in B is 
kept at a constant total energy U, and at constant reservoirs temperatures T±. In this case 
we call Q + , Q- , Qo the work per unit time performed by the forces a+X+ P p a -X- P p a X Ir- 
respectively. 

Let Ce, £+, £- , denote, respectively, the work per unit time performed by the field 

E or by the particles in the thermostats 1Z+,1Z- on the gas in B, or by the gas in B on the 
thermostats 1Z+,1Z-. Then the imposed conservation laws give: 

- Q +C E + £+ + £-= <— ► (U = const) 

2 

-Q++£° =0 ^ ( % ± = N + k B T + ) 

* — ' Im 

l^ n + (3.3) 

2 

-Q + -£°_=0 ^ ( f L = N.k B T.) 

q 67J_ 
— J 

so that one easily finds (by differentiating (3.2) with respect to time and by applying the 
equations of motion, (3.1)) expressions for a_|_,a_,a: 

_ C E + C+ + C- _ Qo _ C± _ Q± 

E £ . 6B Pj/rn ~ E £ . 6B Pj/rn' a± ~ 2N ± k B T ± ~ 2N ± k B T ± 1 ' > 

and the phase space contraction per unit time, or ^- times the entropy creation per unit time, 
is: 

j(x) = (2N± - l)(a+ + a_) + (2N - I) a (3.5) 

in the configuration x. 

The above model also shares the feature that it is time reversible, and the isometric map 
i( p , q ) = (— p , q ) is again such that i S = S~ 1 i, and i 1 = 1. 

§4: The SRB distribution. 

The chaotic hypothesis of §2 allows us to represent the SRB distribution in a simple form, by 
using Markov partitions, [S]. We consider only transitive reversible Anosov systems, although 
many concepts make sense for more general systems with chaotic attractors. 
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The notion of Markov partition is a mathematically precise version of the intuitive idea of 
coarse graining. We just recall here the main properties of Markov partitions. For a discussion 
of the intuitive meaning and the connection with the coarse graining see [G3]. 

1) A parallelogram will be a small set with a boundary consisting of pieces of the stable and 
unstable manifolds of the map S joined together as described below. The smallness has to be 
such that the parts of the manifolds involved look essentially "flat": i.e. the sizes of the sides 
have to be small compared to the smallest radii of curvature of the unstable manifolds W" or 
of the stable manifolds W£ , as x varies in C. 

Therefore let 6 be a length scale small compared to the minimal (among all x) curvature radii 
of the stable and unstable manifolds. Let A" and A s be small (and "small" means of size <C 6) 
connected surface elements on W™ and W£ containing x. We define a parallelogram E in the 
phase space C, to be denoted by A" x A s , with center x and axes A", A s as follows. 

Consider £ £ A u and rj £ A s and suppose that the point z, denoted £ x rj, such that the 
shortest path joining £ with rj formed by a path on the stable manifold W£ joining £ to z and 
by a path on the unstable manifold W™ joining z to rj, is uniquely defined. This will be so if 
6 is small enough and if A", A s are small enough compared to 6 as we assume, (because the 
stable and unstable manifolds are "smooth" and transversal). 

The set E = A" x A s of all the points generated in this way when £,r) vary arbitrarily in 
A", A s is called a parallelogram (or rectangle), provided the boundaries dA u ,dA s of A" and 
A s as subsets of W™ and W£ , respectively, have zero surface area on the manifolds on which 
they lie. The sets d u E = A u x dA s and d s E = dA u x A s will be called the unstable or horizontal 
and stable or vertical sides of the parallelogram E. 

Consider a partition £ = (Ei, . . . , Ej\f) of C into Af rectangles Ej with pairwise disjoint inte- 
riors. We call d u £ = Uf d u Ej and d s £ = Uf d s Ej\ these are called respectively the unstable 
boundary of £ and the stable boundary oi £ , or also the horizontal and vertical boundaries oi £ , 
respectively. 

2) We say that £ is a Markov partition if the transformation S acting on the stable boundary 
of £ maps it into itself (this means Sd s £ C d s £) and if, likewise, the map S 1-1 acting on the 
unstable boundary maps it into itself (S~ 1 d u £ C d u £). 

The actual construction of the SRB distribution then proceeds from the important result of 
the theory of Anosov systems expressed by a remarkable theorem (Sinai, [S]): 

Theorem: Every Anosov system admits a Markov partition £ . 

Comments: 

a) If the reversibility property holds it is clear that i£ is also a a Markov partition. This follows 
from the definition of Markov partition and from the fact that reversibility implies: 

= iw? (4.i) 

b) The definition of a Markov partition also implies that the intersection of two Markov par- 
titions is a Markov partition, hence it is clear that if a transitive Anosov system is reversible 
(i.e. there is an isometry i such that iS = S~ 1 i and i 2 = 1) there are Markov partitions £ that 
are reversible in the sense that £ = i£ , i.e. if Ej £ £ there is j' such that iEj = Ey £ £. 

c) The usefulness of the Markov partitions comes from the possibility that they provide of 
representing the points of C as infinite strings of symbols (in a more useful way than representing 
them, for instance, as the strings of digits that give the value of their coordinates). 

This is simply achieved by associating with x £ C the string j = {jk}'j? = - co , such that 
S k x £ Ej k . The invertibility of the map between x £ C and the compatible or allowed sequences, 
i.e. the sequences j such that the interior of SEj k intersects the interior of Ej k+1 is a well known 
(and easy) consequence of the definition of Markov partition. The correspondence is in fact one 
to one with some obvious exceptions: namely to each sequence j with the above compatibility 
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property corresponds one x; viceversa, if x is not on a boundary of some of the E £ £ nor on 
the image of a boundary under a power of S then x admits only one symbolic representation. 
The points on the boundaries or visiting, in their evolution, the boundaries of course have 
several (but finitely many) symbolic representations. Just in the same way as the decimal 
representation of a number is unique for most numbers: the ones which end with infinitely 
many successive 9 admit two representations. 

The correspondence j< — >x between points x £ C and their history, or symbolic representation, 
j as a compatible sequence will be denoted x = x(j) (symbolic code). 

d) If we define the compatibily matrix, or intersection matrix, C'ij by setting C'ij = 1 if the 
interior of Ej intersects the interior of SEi and C'ij = otherwise, then the assumed transitivity 
implies that there is a iterate q of C such that all elements of C' q are positive (i.e. S q Ej has 
interior intersecting Ek for all pairs j, k, simoultaneosly). 

e) Consider the partition £m = n^ M -5 ,_J '^ obtained by intersecting the images under S k , 
k = —M, . . . , M , of £ . Then £m is still a Markov partition and it is time reversal invariant if £ 
is (see (b) above). Note that the parallelograms of £m can be labeled by the strings of symbols 
j-M , ■ ■ ■ , 3m and they consist of the points x such that S k x £ Ej k for — M < k < M . In other 
words the parallelograms consist of those points x which, in their time evolution, visit at time 
k the parallelogram jj. . 

f) If F(x) is a function on phase space ( observable) then we can regard it as a function F(x( j )) 
on symbolic sequences. An observable F is local if it "depends exponentially little" from the 
history symbols jj. with large k: i.e. if F(x(j)) — F(x(j')) tends to zero exponentially fast 
with the maximum number k such that j and j ' agree on the sites with label h with \h\ < k. 
A simple condition on F guaranteeing its locality is that F is Holder continuous in x. 

We now construct a probability distribution on C by defining it as a probability distribution 
on the space of the compatible strings j and then by interpreting it as a distribution on the 
phase space C. 

(1) For this purpose we first pick a point, that we call the center, %j_ M ,...,j M in each Ej_ My ^j M 
with non empty interior simply by considering the compatible string which is obtained by 
continuing the string j-M, ■ ■ ■ ,3m "to the right" into a string Jm,3m+i, ■ ■ ■ an d to the "left" 
into a string . . . , j-M-i, j-M m a wa Y that the whole string j is compatible (i.e. such that 
there is a point x such that S k x £ Ej k for all k) and, furthermore, the entries of the continuation 
strings depend only on the value of ju and j-M respectively. 

In general given j±M there will be many choices of the continuation strings: which one we 
actually take is irrelevant. We impose however the further constraint that the continuation is 
made in a "time reversible way", i.e. we choose the continuations so that if x is the center of 
Ej then ix is the center of iEj . A further restriction (not necessay in the following but very 
nattural) that one could consider is imposing that the continuation string to the right of jm 
(or to the left of j-M) agree identically after finitely many steps. Note that the existence of 
the continuation strings and the possibility of imposing on them the above restrictions is an 
immediate consequence of the transitivity property of the compatibility matrix C . (2) We then 
define, given r > 0: 

r/2-l 

KA X )= II A -(S j x) (4.2) 

j = -r/2 

where A u (x) is the local expansion coefficient of the surface elements of the unstable manifold 
at x, i.e. it is the jacobian determinant of the transformation S regarded as a map of W™ 
into Wg x . Likewise we define A s (x) and A S T (x) as the corresponding quantities obtained by 
regarding S as a map of the stable manifold W£ to W s Sx . 

(3) Finally we define a distribution /j,m,t on C by "giving" to each set Ej_ My ^j M £ £m a 
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probability proportional to A u 1 T (xj_ Mt ...j M )S T 1 (xj_ M ,...,j M ) where 6(x) = sm-d(x) = 6o(x) is 
the sine of the angle between the stable and the unstable manifolds at x and h T (x) = h(S T l 2 x). 
More precisely we define the distribution /j,m,t so that the integral of a smooth function F is: 

HM,T{dx)F{x) = — J - =rT — (4.3) 

where j is a short hand notation for ]-m, ■ ■ -,3m an d xj = %j_ M ,...,j M is the "center" chosen 
above in Ej £ £m- No relation is assumed here between T and r, although in the applications 
we shall (naturally) take T = r/2, as this simplifies the discussion considerably. 
The distribution /j,m,t is very interesting because it is an approximation (a very good one) 
of the SRB distribution. In fact in [G1,CG1,CG2] the following theorem is shown to be a 
reformulation (convenient although trivially equivalent) of a basic theorem by Sinai: 

Theorem: If (C, S) is a transitive Anosov system the SRB distribution /j, exists and the /j, average 
of a local function (see (f) above) F is: 

n(dx)F(x) = lim / n M , T {dx)F(x) (4.4) 

Furthermore in (4-3) the factor 8~ 1 (xj) could be replaced by 6^(xj) with z any prefixed real 
number (e.g. z = 0). The limits can be interchanged. 

The original statement is that fi exists and it is a Gibbs state with potential log A" 1 (a;): see 
[Bo],[Rl],[R2],[S] for a discussion of this form of the statement. In [R2] the latter statement 
is extended to cover the case in which (C, S) has a global transitive Axiom A attractor: the 
discussion in [G1],[G2] shows that the above theorem extends, unchanged, to such cases. The 
extra factor 6* with z = — 1 was absent in [GC1],[GC2] where z was chosen equal to (an 
admissible alternative choice). 

The possibility of fixing z arbitrarily, in spite of the apparently strong modification it intro- 
duces, is easily seen by examining the proof of (4.4), see [Gl], [G2]. The proof is based on 
the interpretation of (4.3) as a probability distribution on the space of the compatible strings. 
In this interpretation one immediately recognizes that (4.3) corresponds to a finite volume 
Gibbs distribution for a suitable short range hamiltonian defined on the space of compatible 
strings. An extra factor b z T (xj) corresponds to considering the same Gibbs distribution just 
with a different boundary condition: which becomes irrelevant in the limit as r — ► oo because 
one dimensional Gibbs states with short range interactions do not have phase transitions and 
therefore are insensitive to changes in the boundary conditions. Different choices of the center 
points also correspond to different choices of boundary conditions. 

The choice z = — 1 is much better than z = because it leads to simpler formulae and 
arguments: we shall call (4.3) a balanced approximation to the SRB distribution because as 
we shall see it is reminiscent of a probability distribution verifying the detailed balance (which 
however is not verified in our models, except in forcing, i.e. in equilibrium). 

In (4.4) with T ^> r/2 the choice of the point Xj in the parallelograms of £m can be really 
arbitrary, and it does not matter that Xj is really chosen as said above or just anywhere in 
Ej_ Mt ...j M (because the variation of the weigths (4.2) is in this case neglegible, provided M — 
r/2 — ► oo fast enough). 

The extra properties that we need are that £m is reversible (see above), i.e. iEj = Ey £ £m 
(for a suitable j') and that, as a consequence of the reversibility (via (4.1) and the isometric 
nature of the time reversal map i and the validity of y T (%) = —y T {ix) for the underlying 
differential equations): 

A u , r (ia;) = A s \{x), 6 (x) = 6 (ix), 6 T (x) = 6_ T (ix) (4.5) 
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which are identities (see [GC2]); for the definitions of A s ,A s r , see the lines following (4.2). 
Furthermore the volume measure and the expansion and contraction rates are related by: 

A UiT (x)A SiT (x)^\ = e~ Tt °~^ (4.6) 

0—T yX j 

where to is the average time interval between successive timing events and the phase space 
volume contraction for a single transformation is written e~ toa ^ , thus defining a(x) and ~a T (x): 

1 r/2 

<f T (x)= f - < STx ) ( 4 - 7 ) 

r = -r/2 

The (4.6) is is obtained from the relation A u (x)A s (x) = e~ to "^ by evaluating it on the 
points S k x, k = , . . . , ^ and multiplying the results. 

If the time interval t(x) between the timing event i£C and the successive one is very small 
and if its fluctuations can be neglected toghether with those of j(x) see (2.3), (3.5) (within 
the same time interval) then one simply has a(x) = j(x). Note that in all cases with any 
reasonable definition of timing events the time t(x) will tend to zero in the thermodynamic 
limit (as 0(N~ 1 )) ) but also 7 will tend to infinity as O(N). 

More generally there is a simple relation between the function a(x) above and the function 
j(x) which describes the phase space contraction rate in the differential equations giving rise 
to the map S (see (2.3), (3.5)); namely: 

= T / l(S t x)dt (4.8) 
'0 Jo 

But the use of (4.8) is quite clumsy and one can always think that the timing events are chosen, 
artificially, much closer than the natural to = 0(1/ N) and observed at constant time intervals 
so that no difference really exixts between j(x) and a(x). If necessity arises one can always 
use the precise relation (4.8), at the expense of some formal algebraic complications in the 
intermediate step of our coming deductions. 

§5: Applications to the models. Onsager reciprocity and fluctuation dissipation 
theorem. 

In this section we neglect for simplicity of exposition the difference between a(x) and j(x), 
i.e. we suppose that t(x) = to is constant and that 7 is constant on the path traveled in the 
time interval t(x): this simplifies considerably the algebra and the reader should have no trouble 
checking that the proper relation (4.8) could be consistently used leading to no corrections to 
the final results below (because in the end we shall set E_ = 0). 

Relation (4.3) has the form of a statistical average and we shall try to use it in the "same" way 
as in equilibrium statistical mechanics. We shall first study here the two currents Jh, h = 1,2, 
generated by the pair of fields Eh in the diffusion model of §3. 

The two currents are, if ph, i>h are the density and average velocity of the species h : 

J h = PhVh = — ^— = ——— dEh(T (5.1) 

where j £ {h} means that j is a species h particle, and quantities of 0(7V -1 ) have been neglected 
(see (2.3)) and a is j— times the entropy production rate, see (4. 6), (4. 7) and (4.8). Hence if we 
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p ■ 

define T(x), for each configuration x, by Ej ^ = NkT(x) and: 

. de/ 2# J h 1 Ej A tl ~r(*i)^" 1 (*i)c , £;^r(*i) 

where ^(a;) = ^ Er=-r/2 (T ( 5 "' ;C )> see ( 4 - 7 )- If we reca11 ( 4 - 3 ) with | r = M 

we see that can 

be regarded as the SRB average of j^, h = 1,2. 

This expression is similar to the formula derived from the generating function of the Helfand 
moments in [GD2], [GDI]: but it is not the same because in [GD2] the SRB is represented by 
using the notion of (e, r)-separated sets (which are a somewhat more primitive or less concrete 
version of the parallelograms of the Markov partitions). 

We shall also define l u , T ,h,T as: 

A~l(x)8 T (x) - 1 = e T '""W, A SiT (x)8^(x) = e T, '-« (5.3) 
so that (4. 5), (4. 6) imply l UT (x) — l ST (x) = Tto"a T (x). Hence we see that, if dk = dE k '- 
1 Ej ^u^( x i)K 1 ( x i){9hkO : T{xj) + Td k l Ut r(x j )d h a : T (xj)) 



dkjh 



Ej K,l( x j) 6 r 1 {xj) 

1 E; Klixj^^ixi^dkluTixj) E,' K.lixj^-^x^dhWrixj) 



1 T 

-j^{dhk^ T ) - -j-^{{dklu,Tdh^ T ) + {dklu,T){dhV T )) (5.4) 



Here we have interpreted the derivatives with respect to Eh ofa T (x) and l UT (x) by regarding x 
as E_ independent. However this is not quite correct: in fact it is clear that we must consider 
such functions as defined on the attractor, not on the full phase space. The attractor depends 
on £: it can in fact be identified with the unstable manifold Wq of a fixed point O or of a 
periodic orbit O (see [GC2], §4): hence the point xj, which has to be thought as a point on 
the attractor, will change with E_ even though it keeps the same symbolic representation (note 
that the Markov partition changes with E_ although the compatibility matrix does not, by the 
structural E_ stability theorem of Anosov, [AA], at least if E_ is small). 

In taking the derivatives with respect to E_ of l UT (xj) and in defining the current as dE^ri^j) 
there are therefore additional contributions proportional to dE h x. The latter quantity can be 
considered as a function of the symbolic history of the function dE h x(j_) and in [G4] 

it is conjectured that: 

Conjecture: The function dE h x(j) is a local function in the sense of the second theorem in §5 
for all Anosov systems, or axiom A systems, depending smoothly on parameters E_. 

Assuming the validity of the conjecture and using it to perform rigorously an interchange 
of limits one can check, see [G4] for details, that the extra terms in the -derivatives of 
~&T{x{ j )), l UtT (x(j )) at fixed history j just discussed give no contribution to the end result, 
i.e. they do not alter the validity of Onsager's reciprocity or of the fluctuation dissipation 
relation, derived below. Therefore, to avoid formal intricacies, we shall not take into account 
the extra terms and we proceed by ignoring them in (5.4) as well as in the following. The above 
conjecture has a mathematical nature and I do not discuss its proof here: I have not attempted 
to prove it (it seems closely related to Anosov's structural stability theorem, see [AA]). 
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By using the time reversal invariance we see that: 
7^ A,, ^6~ 1 dkl„ r dh7f T 

{Oklu,rOhO- T ) = = (5.5) 

_ J2j {^u^i)^ 1 ^]) dkl u , T {xj) d h a T (xj) +A u 1 T (ixj)6- 1 (ixj)d k l UtT (ix j )d h o : T (ixj j) 
~ 2Z 

where Z denotes the "partition sum" , i.e. the sum in the denominator of (5.4), and the averages 
are with respect to the distribution h t /2,t- 

Recalling that (see (4.5), (4.6)) l UtT (ix) = l S T (x), W T (ix) = —7t t (x) this becomes: 

Ej (K.lixj^r'^XjWklu.Tixj) - K StT {xj)6Z^{xj)d k l StT {xj)) d h W T (xj) 

~z (5 - 6) 

The derivatives at E\ = E2 = can be computed immediately by noting that in such case, 
A UiT (x)A SiT (x) / r ^2) = 1 (see (4.6)). If we use that (4.6) implies l UT — l ST = tofa T then it 
follows, from (5.6), that: 

{{dklu,T d h a T ) - (9fc/ tl)T )(5ftO r r )^ _ = ^ ^(d k a T d h a T ) - (d k a T )(d h a : T )^ (5.7) 

We also see that (since (dhk^r) vanishes in the present case): 

r/2-l r/2-1 

= lim V V ((d k a(S m -)d h a(S m -))-(d k a(S m -))(d h a(S m -))) (5.8) 

m = — r/2 n = — r/2 

where the averages in the r.h.s. are with respect to h t /2,t- 
Hence we see that, apart from a further problem of interchange of limits (see below), (5.8) 
becomes: 

oo 

^'"=212 E (<3^(5 m -)3^(-)>-<3^(-)><3^(-)>) (5-9) 

m = — oo 

where the averages are with respect to the SRB distribution (i.e. to the limit of /U r /2,r)- 
The problem of interchange of limits is easily solved: under our assumption that the system is 
a transitive Anosov system the correlations of smooth observables decay exponentially (because 
they become local observables in the symbolic dynamics interpretation of the evolution, provided 
by the Markov partitions), not only for pL but also for p T /2,T ( m the natural sense in which this 
may be interpreted in a finite r case; e.g. by regarding the interval [—■§■, ■§■] as a circle), and 
uniformly in r. 

The relation (5.9) implies that setting L^k = (dhjk) \ E = o then: 

Lhk = L k h (5.10) 

follows. Note that (5.9) expresses the fluctuation dissipation relation between the transport 
matrix L and the current-current equilibrium correlation. 

In the case of the model in §4 the situation is similar. If ^ . = 2NksT(x) and J q denotes 
the heat — q + = — Q + / L 2 received by the gas from the thermostat 1Z+ per unit time and volume: 

J 2N-l)J2j • H m 1 n 



knT L 2 p 2 /m L 2 

(5.11) 



J_ 2N+ - 1 -Q+ _ _1 

T+ T+ ~ L 2 2N+k B T 2 ~ L 2 



— d ST <T 
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Hence the above argument yields, for the model in §4: 



ST,E = 



ST,E = 



(5.12) 



In general we can consider changing two parameters denoted a, b, thermodynamic forces, which 
control equations of motion of the system. Suppose that the entropy generation per unit time 
a has the form: 



D r 



Qr 



E p]l m 



(5.13) 



where ^ r denotes that the coordinates q are coordinates of a particle belonging to a group 
of N r particles whose phase points are constrained by the r-th constraint that we impose on 
the system (to fix the coordinates that would evolve with a zero Lyapunov exponent, in the 
thermodynamic limit). And let D r be the number of degrees of freedom of the r-th group 
of particles (in 2 space dimensions D r ~ 2N r ); then the above argument can be immediately 
generalized to yield that the flows J a = {d a <r) and Jj = (Sjcr) verify: 



dhJ a 



def 



a,b = 



L21 — daJb 



a,b = 



(5.14) 



which is a general Onsager reciprocity relation between "thermodynamic forces" and "currents" . 
From (5.7) we also see that the matrix <9j J a is positive definite. 

Note that, as mentioned above, in defining d a <r, <9j<7 one has really to think of a as defined on 



the space of the symbolic sequences a 
conceptually different from the "naive" 
affect the end result). 

§6. Remarks. 



this 



<t(x(i)) (so that d a a = j)) + 

(x( j)) although (in the above models it does not 



da 



da dx (Q . 

dx da 



(1) The models in §3, §4 have been considered as undergoing transformations at constant energy 
U . This is not very satisfactory as one also, and mainly, wants to understand also transforma- 
tions in which the internal energy is allowed to change. According to the ideas in §2 this case 
can be treated by imposing that U is constant; but the constant value is fixed on the basis of 
the equation of state of the system. The latter is the relation linking U to the other system 
parameters: 

U = f(E\E 2 ), U = f(E,T-,T+) (6.1) 

in the cases of the models of §3, §4. Here / is determined by the dynamics itself, but its 
computation requires mathematical difficulties that we cannot expect to be able to solve (in 
general). 

Nevertheless we can proceed as in the previous section: in taking the derivatives with respect 
to the parameters there will be extra terms that arise from the fact that the energy surface 
changes as described by (6.1), but the basic symmetry L^k = L^h remains, as one can check 
(by taking also into account that in the equilibrium state obtained when the thermodynamic 
forces, i.e. the variations from the equilibrium values of the parameters, are set to then the 
currents vanish). 

(2) It is quite clear that the discussion of the previous sections can be extended to many other 
models. But it is not clear how far one can really extend the considerations. For instance 
it would be desirable to extend, if possible, the considerations to a microscopic model of a 
macroscopic continuum obeying the macroscopic equations for a fluid or a mixture of fluids 
(possibly with chemical reactions allowed), as defined in [DGM]. 
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(3) Onsager relations are often regarded to be consequences solely of the reversibility of the 
equilibrium dynamics, (see, however, [DJL]). Therefore they must hold also for our models sim- 
ply because they could be derived "as usual", see [DGM], p. 100-101, for a "usual derivation". 

Hence it is worth pointing out that the "usual derivation" rests on several assumptions, none 
of which is needed if the chaotic hypothesis is retained, at least in a dynamics of the type 
considered in the above models. 

(4) The "usual derivation" assumptions are the following; 

(a) linear regression, i.e. the "time behaviour of the state parameters can be described by linear 
equations", see [DGM], p. 36 and p. 100. 

(b) The "Boltzmann postulate" : i.e. the equilibrium distribution is such that the "state param- 
eters" a obey a gaussian large deviation law, see [DGM], p. 91 eq. (46). This means that they 
have an equilibrium probability distribution p( a) proportional to e 2k B G - — = e s (2.)/ k B where 
S is the "entropy" of the state with state parameters a. The entropy is defined by the above 
formula but it is treated as if it had the properties we may expect from a genuine equilibrium 
entropy function (note that this is an assumption). By "large" one means an amount much 
larger than the root mean square values at equilibrium (see [DGM], p. 100), i.e. macroscopic 
(but still very small). One should remark that the gaussian nature is not a consequence of a 
normal distribution assumption as the latter usually concerns small deviation of the order of 
the root mean square values. On the other hand a careful examination of the proof shows that 
all is really needed is that —S( a) is a convex smooth function near a = . Therefore the really 
strong part of the assumption above is that S( a) behaves as an ordinary entropy function (a 
concept that would require some more precise definition). 

(c) The equilibrium evolution is reversible. 

Then it follows that the time evolution of a fluctuation is such that the "state parameters" a 
verify a = LI with X_ = -G a and the symmetry L = L T , see [DGM], p. 101-=-102. 

An initial (distribution of) microscopic configurations, close to the equilibrium state, generates 
a macroscopic state in which fluctuations are possible: so one can consider the free evolution 
of a state in which the initial "state parameters" have an average value off by a from the 
equilibrium value . The evolution will then verify the above "symmetry" relation. 

For instance a = (a\, 02) could be the horizontal center of mass coordinates of the two species 
of particles in model 1 above. 

(5) The connection with "reality" requires further assumptions. Considering our model 1 for 
definiteness, suppose that we act on the system with small forces thus driving an evolution 
of the average values of the "state parameters" t — ► a(t), and creating an entropy per unit 
time (^^+a 2 N 2 E-) ( 

see §2: this makes sense for small fields when the temperature T can 
be identified with the average kinetic energy per particle). Note also that, by (b) above, the 
entropy creation rate is in general S = —G a ■ a = X_ ■ a with X_ = —Ga. 

(d) Then Onsager supposes (see [O]: "As before we shall assume that the average regression of 
fluctuations obey the same laws as the corresponding macroscopic irreversible processes") that 
a fluctuation forced by external forces evolves as if it had occurred spontaneouly; i.e. if X_ is 

v (N 1 E 1 ,N 2 E 2 ) u - 1 1 1 i\t , r N -i • ■ N r E , ■ N r E ■ 

given, X_ = - — - — 7^ — recalling that iVi = N2 = y> ^ ls Q- = T 7" \ 01 1 = T ™> 111 
the notations of the present paper), with Lyi = £21- Note that this is done (and can only be 
correct) up to corrections 0( E_ 2 ). 

Other derivations are more fundamental and are based on kinetic theory (see [C], [DGM]) 
or on the pure microscopic dynamics ([DGM], [D]), but still they rely on various assumptions 
besides time reversibility of the equilibrium dynamics. 
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(6) With our chaotic hypothesis all the above assumptions (a)^-(d) are not necessary, if one 
considers the models in §2, §3, because the final result (i.e. j proportional to L E_ with L 
symmetric and positive definite) has been drawn without further hypotheses (other than the 
mathematical conjecture in §5). Of course assuming the reversibility of the dynamics even 
in nonequilibrium (close to equilibrium) and the chaotic hypothesis is in some sense much 
stronger than the assumptions (a)^-(d) (but note that the reversibility in nonequilibrium , close 
to equilibrium, is a form of the assumption (d)). 

It has, however, the basic advantage of being a conceptually simple general assumption for 
nonequilibrium statistical mechanics, which should furthermore be valid without even the re- 
striction of being close to equilibrium. 

As a final remark one should add that, although reversibility of the nonequilibrium dynamics is 
assumed in this paper, it is quite likely that the ideas and methods can be extended to genuinely 
non reversible models. Attempts at such applications can be already found in [GC1],[GC2],[G4] 
and I hope that they can be generalized to more general physical situations quite beyond the, 
so far special, cases mentioned. 

A much debated question is the extension of the Onsager relations to the really nonequilibrium 
regime (i.e. not close to equilibrium). In this case the very notion of reciprocity becomes ill 
defined: recently the proposal of interpreting the relations as the statement that the flow of 
the state parameters is, in a suitable sense, a "gradient flow" has emerged, [GJL]. The chaotic 
hypothesis in principle applies also in the nonlinear regimes but it is unlikely that it can be 
as powerful a tool to cover deterministic versions of the concrete, exactly treated, stochastic 
models of [GJL]. 



Appendix Al:The Gauss minimal constraint principle. 

Let ip( i , x ) = be a constraint and let R(i,x) be the constraint reaction and F^(i, x) 
the active force. 

Consider all the possible accelerations a compatible with the constraints and a given ini- 
tial state i , x . Then R is ideal or verifies the principle of minimal constraint if the actual 
accelerations a i = -^-{ F_ i + R { ) minimize the effort: 

JV N 

■ 1 m i ■ 1 

z=l z=l 

for all possible variations 6 a i compatible with the constraint <p. 
Since all possible accelerations following x_ , x_ are such that X^fci ®& ■£.) ' ^ £U = we 

can write: 



F_ i — rrii a i — adi f(x_, x) = (Al.2) 



with a such that -^<p(x_, x_) = 0, i.e. : 

a = — j 

which is the analytical expression of the Gauss' principle. 
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